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ABSTRACT 

We report the strength of seed magnetic flux of accretion disk surrounding 
the PopIII stars. The magnetic field in accretion disk might play an important 
role in the transport of angular momentum because of the turbulence induced by 
Magneto- Rotational Instability (MRI). On the other hand, since the primordial 
star-forming clouds contain no heavy elements and grains, they experience much 
different thermal history and different dissipation history of the magnetic field in 
the course of their gravitational contraction, from those in the present-day star- 
forming molecular clouds. In order to assess the magnetic field strength in the 
accretion disk of PopIII stars, we calculate the thermal history of the primordial 
collapsing clouds, and investigate the coupling of magnetic field with primordial 
gas. As a result, we find that the magnetic field strongly couple with primordial 
gas cloud throughout the collapse, i.e. the magnetic field are frozen to the gas as 
far as initial field strength satisfies B < 10~ 5 (?t,h/10 3 cm" 3 ) 0,55 G. 

Subject headings: accretion, accretion disks — diffusion — early universe — stars: 
formation — stars: magnetic fields 



1. INTRODUCTION 



PopIII stars are considered to have very important impacts on the thermal history and 
the chemical evolution of the universe. Recent observations of the polarization of Cosmic 
Microwave Background (CMB) photons by WMAP (Wilkinson Microwave Anisotropy Probe) 
(Kogut et al. 2003) revealed that the optical depth of the universe by Thomson scattering is 
fairly large (r = 0.17 ± 0.04). This result tells that reionization epoch is earlier {z — 17 ± 5) 
than expected from the observations on Gunn- Peterson trough (Becker et al. 2001; Fan 
et al. 2002; White et al. 2003). On the other hand, the results of theoretical calculations 
of structure formation (Sokasian et al. 2003) indicate that such early reionization seems 
to be difficult solely by the Pop II stars, but it might be possible if the ionizing photons 
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from PopIII stars with top heavy initial mass function (IMF) is taken into account. Thus, 
the mass (or IMF) of PopIII star have great significance on the reionization epoch of the 
universe. In addition, some amount of metals are found at high redshift intergalactic matter 
by the observation of high redshift QSO absorption line systems (Songaila 2001; Ma 2002; 
Vladilo 2002). This means significant fraction of baryons are already processed in the stars 
by z = 5 (Songaila 2001). PopIII stars also should be responsible for such metal pollution 
of the early universe, and the abundance pattern of heavy elements depends on the mass of 
PopIII stars. Thus, typical mass (or IMF) of the PopIII stars is a key quantity also for the 
chemical evolution of the universe. 

In order to assess the typical mass of PopIII stars, typical scale of the prestellar core 
(or fragments of primordial gas) should be evaluated as a first step. Several authors have 
studied on this issue by simple one-zone approach before middle of '90s (Matsuda, Sato & 
Takeda 1969; Hutchins 1976; Carlberg 1981; Palla, Salpeter & Stahler 1983; Susa, Uehara & 
Nishi 1996; Uehara et al. 1996; Puy & Signore 1996). Recently, multi-dimensional numerical 
simulations of fragmentation of primordial gas have been performed intensively by several 
authors(Nakamura & Umemura 1999; Bromm, Coppi & Larson 1999; Abel, Bryan & Norman 
2000; Nakamura & Umemura 2001; Bromm, Coppi & Larson 2002), and they find that the 
mass of PopIII prestellar core is quite massive (~ 10 3_4 M Q ). 

Further evolution of prestellar core was investigated by Omukai & Nishi (1998), and it 
is found that the collapse proceeds in a run-away fashion and converges to Larson-Penston 
similarity solution (Larson 1969; Penston 1969) with 7 ~ 1.1. They also find the mass 
accretion rate is also very large compared to the present-day forming stars, although spherical 
symmetry is assumed in their numerical calculation. In reality, however, prestellar cores 
have some amount of angular momentum, which prevent the mass from accreting onto the 
protostars. Consequently, accretion disks surrounding the protostars are expected, and some 
mechanisms that transport their angular momentum are required in order to enable the mass 
accretion onto protostars. 

There are a few possibilities of the angular momentum transport such as gravitational 
torque by the nonaxisymmetric structures in the accretion disk, the interaction among the 
fragments of the disk (Stone et al. 2000; Bodenheimer et al. 2000), and the turbulent vis- 
cosity triggered by Magneto-Rotational Instability (MRI) (Hawley & Balbus 1992; Sano, 
Inutsuka & Miyama 1998; Sano & Inutsuka 2001). All of these mechanisms are regarded to 
be important for present-day star formation. Thus, it is worth to investigate these processes 
for primordial case. Among these possibilities, we concentrate on the last one, the MRI 
induced turbulence. 

In order to activate MRI, the initial magnetic field strength in the accretion disk should 
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be larger than a critical value. Otherwise the magnetic field is dissipated before the field is 
amplified by MRI (Tan k McKee 2004; Tan & Blackman 2004). Thus, it is important to 
assess the "initial" magnetic field strength in the accretion disk. As was discussed by Nakano 
& Umebayashi (1986) for preset-day case, magnetic field could be dissipated while prestellar 
core collapses. However, the dissipation processes strongly depend on the components of the 
gas as well as the the temperature evolution in the course of the collapse. On the other hand, 
the temperature of the collapsing primordial gas is much higher than that of the present-day 
case (Omukai 2000), which might bring about different dissipation history of the magnetic 
field. 

In this paper, we investigate the dissipation of magnetic field in collapsing primordial 
gas cloud. Consequently, we obtain the initial magnetic field strength in the accretion disk of 
PopIII stars. In the next section, formulations are described. In §3, results of our calculations 
are shown, and the importance of MRI in PopIII star formation is discussed in §4. Final 
section is devoted to summary. 



2. CALCULATIONS 

We consider the collapse of spherical clouds without metals and dusts, but with slight 
magnetic fields. Then, we calculate the time evolution of the central density, chemical com- 
position and dissipation of magnetic flux from the gravitationally collapsing core. The initial 
strength of magnetic fields when primordial clouds begin to contract is given as unknown 
parameter. 



2.1. Dynamics 

We assume that the dynamics is described by the free-fall relation, 

dp_P_ m 

dt T ff ' 1 } 

where p is the density of the collapsing core and Tg = (37r/32Gp) 1//2 is the free-fall time. 
In fact, as described in Omukai (2000), thermal evolution of the core in ID hydrodynamic 
simulation is well described by such one zone model. On the other hand, we also have to 
evaluate the dissipation of magnetic field in the accretion phase, and it is not the same as 
the case of collapsing core. In this case, the temperature of the accretion flow would rise 
faster than those in the core, because of shock heating. Thus, the magnetic field is expected 
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to be frozen to the gas also in the accretion flow, if the frozen condition is satisfied in the 
collapsing core. 

It is also worth noting that this free-fall approximation is also marginally valid for the 
case of quasi-hydrostatic core found in the simulation by Abel, Bryan & Norman (2000, 
2002). In fact, panel C of Fig.3 in Abel, Bryan & Norman (2000) tells that the cooling time 
is as short as the free-fall time. Thus, the system evolves with time scale r ff ~ t coo i. 

In order to describe thermal processes, we approximate the equation of state to be 
polytrope with index 7 = 1.1, 

V = Kp\ (2) 

where K is the constant coefficient. The reason to adopt the polytrope of 7 = 1.1 is that 
in past work the thermal evolution of primordial collapsing core is investigated in detail by 
Omukai (2000), where the temperature grows with 7 ~ 1.1 in considering heating/cooling 
processes. 

Since we are interested in the collapsing gas cloud, the magnetic energy density needs 
to be less than the gravitational energy density of the cloud core. The critical field strength 
B CT is defined by the equation B^/AtvR ~ pGM/R 2 . Since the radius of the core should be 
on the order of Jean's length, B cr is given by the following equation, 

B cr ~ (3G) 1/2 Mj/i23 

= (4/3) 1 /V/ 2 c sP 1 / 2 oc p - 55 , (3) 

where Rj, Mj, and c s are the Jean's length, the Jean's mass, and the sound velocity, respec- 
tively. 

2.2. Dissipation of Magnetic Fields 

There are two processes for the dissipation of magnetic fields. 1 The first one is Ohmic 
loss, and the other is ambipolar diffusion. We assess the degree of dissipation defined by 
the ratio of the drift velocity v-q x of the field lines from the gas to the free-fall velocity, as 
was introduced in Nakano & Umebayashi (1986). In Nakano & Umebayashi (1986), the drift 
velocity involves both the Ohmic dissipation and the ambipolar diffusion processes. There 
are two important quantities which characterize these diffusion processes. They are r v and 



1 Another process is the diffusion of unidirectional magnetic fields by turbulent flows (Kim & Diamond 
2002), which is beyond the scope of this paper, since we investigate the frozen- in condition of the directional 
fields. 
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u u which denote the viscous damping time of the relative velocity of charged particle v to the 
neutral particles, and the cyclotron frequency of the charged particle z/, respectively. Then, 
r v is expressed as 

where is the reduced mass, n u , n n , and p v are, the mean number density for the charged 
particle u, the neutral particle n, and the mass density of charged particle u, respectively 
The averaged momentum-transfer rate coefficient for a particle v colliding with a neutral 
particle is expressed by (av) un . 

The momentum-transfer cross section of an ion with a neutral particle is given by the 
Langevin rate coefficient (Osterbrock 1961). For an electron, it is found experimentally by 
Hayashi (1981) that the cross sections at low energies are much smaller than the Langevin 
rate coefficient and are nearly equal to a geometrical cross section. We use the empirical 
formulae for the momentum-transfer rate coefficients (Kamaya & Nishi 2000; Sano, Miyama, 
Umebayashi, & Nakano 2000). 

According to Nakano & Umebayashi (1986), the drift velocity can be given by 

VBx = ^r-{j x B) x , (5) 
A c 



where 



A = A\ + A\, (6) 



1 V V 



(8) 



B is the mean magnetic field in the primordial cloud, the suffix x means x direction com- 
ponent in local Cartesian coordinates where the z direction is taken as the direction of B. 
In calculating the drift velocity vb x according to equation (5) we replace (l/c)(j x B) x with 
the mean magnetic force B 2 /47rR, where B is the mean field strength in the cloud, R is the 
radius of the cloud. 

The magnetic field is mainly dissipated by the Ohmic loss when \t v uj v \ < 1 for main 
charged particles. In such a situation we have the approximate expression from equations 
(5), (6), (7) and (8) that 
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where a c the electrical conductivity, 



q u and m v are the electrical charge and the mass for a charged particle u, respectively. Thus 
the drift velocity is independent of B. On the other hand, when \t v uj v \ > 1, the ambipolar 
diffusion is a main process of dissipation. In such a case we obtain 

n B 2 

t>R T ~ , (11) 

where the suffix i expresses the ion particles. Note that the drift velocity is proportional to 
B 2 in this case. 

Tan & Blackman (2004) have evaluated the drift velocity of magnetic field in the quasi- 
hydrostatic core found in Abel, Bryan & Norman (2000, 2002), when its density is a certain 
value. However, this system is highly non-equilibrium, and the ionization degree gets smaller 
and smaller as the collapse proceeds. Thus, its is not trivial at all whether the frozen 
condition is satisfied or not at much higher density. Therefore, we have to perform detailed 
non-equilibrium calculations of chemical reaction network, in order to obtain the correct 
ionization degree in the collapsing gas. 



2.3. Chemistry 

In order to investigate the evolution of ionized fraction during the collapse in detail, we 
solve non-equilibrium chemical reaction network of primordial gas that involves not only H 
element, but also D, He, and Li. Furthermore, we introduce following 24 species: e~, H + , 
H, H", H 2 , H+, H+, D, D+ D" HD, HD+, H 2 D+, He, He+, He++, HeH+, Li, Li+, Li++, 
Li 3+ , Li - , LiH, and LiH + . We employ the latest reaction rate coefficients appearing in the 
following papers, Galli & Palla (1998), Omukai (2000), Standi, Lepp & Dalgarno (1998), 
Flower (2002) and Lepp, Standi & Dalgarno (2002). As for the radiative recombination, we 
use the rate coefficients based on Spitzer (1978). 



2.4. Initial Conditions 

We employ two different initial conditions A and B. For model A, we set nn.ini = 
10 3 cm~ 3 and T ini = 250 K for initial conditions of prestellar core. We use the fraction of 
the chemical composition in the early universe given by Galli & Palla (1998) as the initial 
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condition for the cloud. This model is on the evolutionary track (Z = 0, fiducial model) 
in Omukai (2000). On the other hand, as was shown by Palla, Salpeter & Stahler (1983) 
and Omukai (2000), the paths in nn — T plane rapidly converge to an identical trajectory 
from various initial conditions. Remark the almost hydrostatic prestellar core found in the 
numerical simulations in Abel, Bryan & Norman (2000) is also very close to the above 
trajectory. Thus, this choice of initial condition is quite natural for the formation of PopIII 
stars. 

On the other hand, there is a case which do not converges rapidly to the trajectory, 
which was pointed out by Uehara & Inutsuka (2000) and Nakamura & Umemura (2002). In 
this case, the gas is initially heated up to 10 4 K during the collapse of rather massive host 
galaxy, and H 2 is formed efficiently utilizing the residual electrons as catalysts (e.g. Susa 
et al. 1998). Afterwards, HD is formed from H 2 , and the gas is cooled down to ~100 K. 
Furthermore, the gas temperature is kept at relatively low level (~ 100 K) in the course 
of the collapse. We employ this case as model B, in which we set the initial condition 
(^H,mi = 10 3 cm~ 3 ,T ini = 100 K), and the initial fractional abundances are taken from the 
results of Nakamura & Umemura (2002). 

The initial magnetic field strength for the primordial star-forming environment has been 
investigated by many authors (e.g. Pudritz & Silk 1989; Kulsrud et al. 1997; Widrow 2002; 
Langer, Puget & Aghanim 2003). However, the initial field generation is still a controversial 
question. Thus, the initial magnetic field strength is given as the model parameter in our 
calculations. 

3. RESULTS 

3.1. Ionization degree 

Figure 1 shows the fractional abundance of various species, e, H + , H, H 2 , Hjj~, Li, and 
Li + for model A. At low density, n H < 10 11 cm -3 , the ionization degree decreases as density 
increases because of the recombination of the electrons with the protons. However, since the 
Li + recombination rate is smaller than that of H + , the reduction of the ionization degree is 
almost quenched at n H ~ 10 12 cm" 3 . Moreover, the ionization degree increases at the higher 
density because the collisional ionization is activated by high temperature. Note that the 
ionization degree is very small, but there is a minimal value ~ 10~ 12 . Figure 2 shows the 
results for model B. In this case, the fraction of electrons evolve in a similar way. There is 
a minimal value ~ 10~ 12 , at which the electrons are provided by Li. 

In order to clarify the importance of Li, we perform control runs without Li, Li + , Li ++ , 
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Li 3+ , Li - , LiH, and LiH + , for both initial conditions A and B. The results of these runs 
are shown in Figures 3 and 4. In the case of model A, even without Li, ionized fraction do 
not gets smaller than ~ 10 _12 (Figure 3). Therefore, presence of Li is not so important in 
this case. However, for the model B, the result is much different from the other case. In 
the presence of Li, electron fraction is kept around 10 -12 even at wh > 10 11 cm~ 3 . But in 
the absence of Li, electron fraction gets much smaller than 10~ 12 (Figure 4). Thus, Li is 
quite important in the low temperature model. We will come back to this point in the next 
subsection. 

We also emphasize the significance of non-equilibrium treatment of chemical reactions. 
Figure 5 shows the fractions of various chemical species after the time integration assuming 
fixed density and temperature over 10 3 Gyrs (a hundred times larger than the age of the 
universe). As a result, for n H > 10 13 cm -3 , the fractions of chemical species converge 
to equilibrium values, although they are still not in equilibrium for lower densities. It is 
clear that results from non-equilibrium calculation is very different from those in chemical 
equilibrium. Especially, the non-equilibrium fraction of electron is much larger than the 
equilibrium value. The time scales of equilibration are also shown in figure 6, and they are 
always much longer than the local free-fall time. Thus, non-equilibrium calculations are 
indispensable to evaluate the ionization degree of collapsing primordial gas. 



3.2. Drift velocity of magnetic field 

We calculate the drift velocity of magnetic field using the results of the time dependent 
ionization degree in Figures 1-4. The results are shown in Figures 7-10 . The dotted curves 
in Figures 7-10 show the contours of the logarithmic value of the drift velocity normalized 
by the free-fall velocity, vb x /ue on n^-B plane. The solid curve represents the locus along 
which v Bx = u s is satisfied. Here us denotes the free-fall velocity. 

The condition |t c cj c | = 1 holds for electrons at the field strength 

B = cfi cn n n (av) cn ^ 
e 

The field strength B e is also shown by the dashed line. When the magnetic field is stronger 
than B e , i.e. |r c c<j c | > 1, the ambipolar diffusion is the dominant process for the field 
dissipation. On the other hand, when the magnetic field is weaker than B c , i.e. |t c cj c | < 1, 
magnetic flux is lost due to the Ohmic dissipation. 

The dot-dashed line represents the critical field B CT . As we mentioned previously, since 
we are interested in the collapsing cloud, our calculations are appropriate in the region 
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where the field strength satisfies B < B cr . It is clear from Figures 7 and 8 that the frozen-in 
condition v-q x < is almost alway satisfied as far as B is less than the critical field strength 
B CT . Therefore, the magnetic field is always frozen to the gas as long as we consider the 
collapsing clouds for both of the initial conditions A and B. 

Since the field is frozen to the cloud, the magnetic flux is conserved during the contrac- 
tion, i.e. 

-2/3 



B 



in- 5 ( B[n[ \ ( nH ' ini V 
VlO- 11 GJ V10 3 cm" 3 / 



x " H - G. (13) 



. 10 12 cm" 

Following the formula of disk radius r<j (equation 17) in Tan & McKee (2004), we can assess 
the density when the disk is formed: 



„ d . 2 . 9xU >*(|f)~ 6 (^)~* m - (14) 

Here M denote the total mass within the radius of the disk, f Kcp is the ratio of circular 
velocity and Keplarian velocity at the sonic point of the accretion flow. Remark that is 
not the density of the accretion disk. It is evaluated by the simple formula 3M/ (47rr^/im p ), 
that is the averaged density within the radius r&. Thus, this density could be compared 
with the density of the core calculated in our model. Combining equations (13) and (14), 
we obtain the magnetic field strength in the accretion disk: 



B = 2.0 x 10~ 5 



Bmi \ / n H ,ini \" 2/3 



10" 11 G VlO 3 cm" 3 



( fa cp \ 4 ( M 



-40/21 

HifJ feJ & (15) 

We also evaluated the diffusion velocities in the runs without Li. Figures 9 and 10 
represents the results for models A and B. As expected from the results found in the previous 
subsection, the diffusion velocity in model A is not so different from the results with Li, 
however, in model B, magnetic field is lost by the Ohmic loss at n H > 10 16 cm~ 3 . Of course, 
this is not the realistic calculation, but we can learn the significance of Li from this results. 



4. Discussion 



The condition for which the MRI can grow at the PopIII accretion disk is investigated by 
Tan k, Blackman (2004). This condition requires that the MRI growth timescale is shorter 
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than the diffusion timescale. Thus, there is the minimum field strength in the disk for the 
MRI to be driven. According to Tan & Blackman (2004), this minimum field strength in 
the disk becomes 



where is the stellar mass, r is the radius from the cloud center, p disk denotes the density 
of the accretion disk and In A is the Coulomb logarithm. Characteristic values of these 
parameters are taken from Figure 3 in Tan & Blackman (2004). 

Hence, it is concluded that if the initial field strength in prestellar core is at least 
> 10~ 10 G at rin = 10 3 cm -3 , the transport of angular momentum could be driven by the 
turbulence due to the MRI in the accretion disk. 

In fact, there are the models for the generation of initial seed magnetic field in the 
universe (Pudritz & Silk 1989; Kulsrud et al. 1997; Widrow 2002; Langer, Puget & Aghanim 
2003). Among these models, Langer, Puget & Aghanim (2003) propose the generation 
mechanism of magnetic field based on the radiation force around very luminous objects such 
as QSOs. According to their study, the magnetic field ~ 10~ n — 10~ 12 G is generated in the 
inter galactic matter. Consequently, field strength is amplified to ~ 10~ 7 — 10~ 8 G when the 
clouds collapses to nn = 10 3 cm -3 , which is the initial condition of our analysis. Thus, in this 
case, the MRI can be driven and it could be the possible mechanism of angular momentum 
transport. On the other hand, most of the other seed field generation mechanism predict 
B < 10~ 19 G, which is too small to drive MRI. Thus, MRI may not be important for the 
formation of very first stars, since the number of sources which provide anisotropic radiation 
field should be very small at the epoch of very first stars. 



In this paper, the dissipation of the magnetic field in the collapsing primordial gas cloud 
is investigated using a simple analysis. As a result, we find that that magnetic field is frozen 
to the gas as far as the initial field strength satisfies B < 10~ 5 (n H /10 3 cm" 3 ) - 55 G. This 
condition holds when the magnetic field does not affect the dynamics of the gravitational 
contraction. It is also found that seed magnetic field induced by radiation force is strong 
enough to activate MRI in the accretion disk surrounding PopIII stars. The MRI induced 



B 




(16) 



5. 



SUMMARY 
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turbulence might play an important role in transporting the angular momentum in PopIII 
accretion disk. 
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Fig. 1. — The evolution of fractional abundance of main species, e, H + , H, H 2 , H3 , Li, and 
Li + , are plotted for model A. Vertical axis denotes the fractional abundance of the above 
each species and horizontal axis is the density of the cloud. 
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Fig. 2. — Same as Figure 1 except that the initial temperature is 100K (model B). 
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Fig. 3. — Same as Figure 1 except that this model is calculated without Li. 
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Fig. 4. — Same as Figure 2 except that this model is calculated without Li. 
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Fig. 5. — The fractions of various chemical species after the time integration assuming fixed 
density and temperature over 10 3 Gyrs are superimposed on Figure 1. Curves with markers 
represent the fractions after long time integration, and the unmarked curves represent the 
non-equilibrium fractions. 



-19- 




Fig. 6. — Time scale of equilibration is shown for several species. 
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Fig. 7. — The drift velocity vb x of the magnetic field as a function of the density nu and the 
field strength B for a collapsing cloud for model A. The dotted curves represent the contours 
of constant v-q x normalized by the free-fall velocity ug and its logarithmic values are labelled 
on those curves. The solid curve shows the locus along which v& x — ug is satisfied. The 
dashed line represents the locus \t v uj v \ = 1 for electrons. The dot-dashed line represents the 
critical field B cr given by equation (3). 
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Fig. 8. — Same as Figure 7 except that the initial temperature is 100K (model B). 
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Fig. 9. — Same as Figure 7 except that this model is calculated without Li. 
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Fig. 10. — Same as Figure 8 except that this model is calculated without Li. 



